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ABSTRACT 

A newly developed anisotropic K-e model is applied to calculate three axisymmetric 
diffuser flows with or without separation. The new model uses a quadratic stress-strain 
relation and satisfies the realizability conditions, i.e., it ensures both the positivity of 
the turbulent normal stresses and the Schwarz’ inequality between any fluctuating veloc- 
ities. Calculations are carried out with a finite-volume method. A second-order accurate, 
bounded convection scheme and sufficiently fine grids are used to ensure numerical cred- 
ibility of the solutions. The standard K-e model is also used in order to highlight the 
performance of the new model. Comparison with the experimental data shows that the 
anisotropic K-e model performs consistently better than does the standard K-e model in 
all of the three test cases. 


1. INTRODUCTION 

Diffusers are widely used in aeropropulsion systems such as engine inlet, nozzle and 
combustion chamber. The maneuverability of modern aircraft at a wide range of speeds 
requires to control engine face flow distortion. Flow recirculation in combustors is critical 
to stabilize flames and to increase the fuel mixing efficiency. Therefore, the understanding 
of diffuser flows is of paramount importance to the design of such devices. Until recently, 
calculations of diffuser flows in the aeropropulsion community have been almost done with 
either mixing-length models such as the Baldwin-Lomax model or two-equation isotropic 
eddy viscosity models. The computational results are often less than satisfactory for the 
purpose of design, and this is mainly attributed to turbulence modeling. The limitations of 
the mixing-length models are quite apparent, and it is questionable whether improvements 
in predicting general diffuser flows can be made using this kind of models. Two-equation 
isotropic eddy viscosity models also have some undesirable features such as isotropy and 
poor response to both adverse pressure gradient and streamline curvature effects. 

Recently, Shih and Lumley (1993) have developed a sixth-order stress-strain relation 
under the assumption that the turbulent stress is a function of the mean deformation 
tensor, the velocity and length scales of turbulence characterized by the turbulent ki- 
netic energy and its dissipation rate. This is the most general constitutive relation for 
the Reynolds stress within the framework of algebraic turbulence modeling, with the lin- 
ear stress-strain relation in the Boussinesq’s eddy-viscosity concept being its first-order 
approximation. The general relation contains eleven coefficients to be determined, an 
inconvenience for practical applications. Therefore, Shih et al. (1994) have developed 
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a new anisotropic K-e model which uses a truncated quadratic stress-strain relation in 
conjunction with the two modeled transport equations for the turbulent kinetic energy 
and its dissipation rate. The new model has the following desirable features: (a) it is 
fully realizable, i.e., it ensures both the positivity of the turbulent normal stresses and 
the Schwarz’ inequality for turbulent shear stresses; (b) it accounts for the effect of the 
mean deformation rate by which the eddy-viscosity will be maintained at an adequate 
level to mimi c complex flow structures; (c) it shows the proper lack of a rotation effect 
on the isotropic turbulence, satisfying the rapid distortion theory; and (d) it is easier to 
use, as compared to other formulations. Simplicity is of great value for practical engi- 
neering applications. The model were calibrated using both well-studied basic flows (e.g. 
homogenous shear flows and the surface flows in the inertial sublayer) and complicated 
backward-facing step flows. 

In this work, the new anisotropic K-e model is applied to calculate three axisymmet- 
ric diffuser flows, two attached and one separated, all of which undergo strong adverse 
pressure gradients. The standard K-e model (Launder and Spalding, 1974) is also used, 
since it is the most popular turbulence model used today in calculations of complicated 
flows and can also be used to highlight the performance of the new model. The calcula- 
tions are carried out with a finite-volume method. The numerical accuracy of solutions 
is ensured by using a second-order accurate, bounded convection scheme and sufficiently 
fine grids. The performances of the models are examined through extensive comparisons 
with experimental data. 

2. CALCULATION PROCEDURE 

Incompressible, steady-state, turbulent flows are governed by the following Reynolds 
averaged continuity and Navier-Stokes equations: 

= o (1) 

(pUjU t — pU x>3 — t, 3 )j = —p t (2) 

where U t is the mean velocity component in indirection, p is the pressure, p and p are the 
fluid molecular viscosity and density, and U x>3 is the derivative of U t with respect to the 
co-ordinate x 3 , respectively. The Reynolds stress r tJ (= —pu x u 3 ) in Eq.( 2) is calculated 
by using the following two turbulence models: 

1) Standard K-e model 
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C„ = 0.09, C u = 1.44, C 2t = 1.92, a K = 1 , <r„ = 1.3 
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2) Anisotropic 

K-e model 
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The K and e in the anisotropic K-e model are calculated with the same equations as in 
the standard K-e model. 

The above two models are of high Reynolds number form, therefore, they cannot be 
integrated down to the wall. In this work, the standard wall function approach (Launder 
and Spalding, 1974) is used to bridge the near-wall region. 

Numerical calculations are carried out by using a conservative finite- volume method 
designed for calculating two-dimensional planar or axisymmetric, incompressible, elliptic 
flows with complex boundaries. The method uses a non-staggered variable arrangement, 
namely, all the dependent variables are stored at the geometric center of each control vol- 
ume. The momentum interpolation of Rhie and Chow (1983) is used to avoid checkerboard 
oscillations usually associated with non-staggered grids. The velocity-pressure coupling 
is achieved via the SIMPLEC algorithm (Van Doormal and Raithby, 1984). Numerical 
implementation of the anisotropic K-e model is straightforward; the linear part in Eq.( 8) 
is treated in the same manner as in the standard K-e model and the quadratic part is 
included in the source terms. To ensure numerical accuracy and stability, the convec- 
tion terms of all the transport equations axe differenced by the hybrid linear/parabolic 
approximation (HLPA) of second-order accuracy (Zhu, 1991a), and all the other terms 
by the conventional central differencing scheme. The resulting set of algebraic difference 
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equations is solved with the strongly implicit procedure of Stone (1968). The details of 
the present numerical procedure are given in Rodi et al. (1989) and Zhu (1991b). 


3. RESULTS 

All calculations were performed on the Cray YMP computer. Grid sensitivity tests 
were conducted using different grids ranging from 51x41 to 161x100, and the results 
to be presented in the following are all grid-independent. The length of computational 
domain was taken as 20Ro with Ro being the radius of the diffuser at entrance. 

3.1 Azad and Kassab’s case 

The experimental apparatus of Azad and Kassab (1989) was a long pipe followed by 
the diffuser with an 8° included angle. The Reynolds number based on the diameter of 
the pipe and the bulk mean velocity in the pipe was 1.15 x 10 5 and the fully developed 
flow conditions were at the entry of the diffuser. The calculations were made using a 
83x53 grid. The first grid point away from the duct wall had the value of y + varying 
between 36 ~ 67. At the inlet of the diffuser, the fully developed flow conditions were 
specified which were obtained from a separate calculation of the flow in a long pipe. The 
iterations and CPU-time in seconds required to get the maximum normalized residue of 
all the dependent variables below 10 -6 were 546 and 74 for the anisotropic K-e (AKE) 
model, 499 and 64 for the standard K-e (SKE) model. 

Figures 1 and 2 show the axial variations of calculated and measured wall static pressure 
coefficient and centerline velocity. It is seen that the results of the AKE model are in 
good agreement with the experimental data along almost the entire length of the diffuser, 
while the SKE model overpredicts the pressure coefficient and underpredicts the centerline 
velocity decay in the downstream region of the diffuser. Figures 3 and 4 show the profiles of 
turbulent kinetic energy and shear stress at three axial positions; here both the turbulent 
kinetic energy and shear stress are normilized by the bulk mean velocity in the pipe. 
Although the agreement with the experimental data is less satisfactory in figures 3 and 
4 for the turbulent quantities than in figures 1 and 2 for the mean flow quantities, the 
AKE model is still seen to performs better than does the SKE model in predicting the 
turbulent quantities, overall. 

3.2 Fraser’s case 

The diffuser in the experiment of Fraser (1958) had a 10° included angle, and the 
Reynolds number based on the inlet diameter and the free stream velocity was 5.08 x 10 5 . 
This flow was also a test case for the 1968 AFOSR-IFP-Stanford Conference (Coles and 
Hirst, 1968). Unlike the Azad and Kassab’s case, the Fraser’s flow was of boundary layer 
type. The computational domain was discretized by a 145x52 grid. No turbulent data 
was available and it has been found that inlet turbulent boundary conditions had a big 
influence on the computed results. In this study, these boundary conditions were obtained 
from the outlet profiles of a pipe flow, calculated separately and by try and error so that 
the flow velocity at the pipe outlet matched the experimental data at the inlet of the 
diffuser. 

Figures 5 and 6 show the friction coefficient and the static pressure coefficient along 
the duct wall. Figures 7 to 9 show three boundary layer quantities, i.e., the displacement 
thickness, momentum thickness and shape factor. Figure 10 shows the centerline velocity 
decay and figure 11 shows the axial mean velocity profiles at three axial positions. It can 
be seen from all these figures that the results of the AKE model are consistently better 
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than those of the SKE model, especially at downstream locations. 

3.3 Binder and Kian’s case 

The flow configuration of the confined jet, measured by Binder and Kian (1983), is 
sketched in figure 12. The diffuser has a half angle 9 = 2.5° and the inlet diameter 
D 0 = 16cm. At the entrance, two uniform flows, a jet of larger velocity Uj from the nozzle 
of diameter d 0 = 1.6cm and an ambient stream of smaller velocity U a , are discharged into 
the diffuser. The inlet flow conditions can be characterized by the Craya-Curtet number 
C t and the experiment showed that recirculation occurs when C t < 1.10. Several C t values 
were considered in the experiment, but in the present study, only the flow at C t = 0.59 
was calculated, which corresponded to Uj = 40cm/s and U a = 2.33cm/s. Under this 
condition, the flow had a strong separation and the previous calculations gave the worst 
agreement with the measurement (Zhu, 1986). A 86x50 grid was used in this calculation. 

Figure 13 shows the variation of pressure coefficients along the duct wall. Here, C p is 
defined by 

Ap-pUl/2 

m 2 

and A p is the pressure difference between the location x and the entrance. The pressure 
variation is governed by the jet entrainment, the contraction and expansion of the flow 
caused by recirculation bubble as well as the geometry of the duct. It can be shown that 
the entrainment and the divergence of the duct can only produce a maximum pressure 
difference equal to pU%/2, while the pressure rise in the recirculation zone depends on 
the width of the recirculation bubble and can be much larger than pUH 2. Regarding 
the comparison between predictions and experiments, it is seen that both models give 
good agreement with the experiment in the upstream region. However, in this region, 
the pressure gradient is mainly determined by the ambient potential flow on which the 
turbulence models have little effect. In the downstream region where recirculation occurs, 
the AKE result is still in good agreement with the experiment, while the SKE model 
predicts a premature pressure rise. 

Figures 14 and 15 show the profiles of axial mean velocity and turbulent shear stress 
at three axial positions; here U m and U 0 are the sectional average velocity and centerline 
velocity, respectively. For the velocity profiles shown in figure 14, the comparison between 
the predictions and experiment is generally good, but both models fail to capture some 
detailed flow features. At the location x)D 0 = 1.25, the experimental profile shows 
no constant ambient velocity portion while the calculated profiles by both models still 
have a vertical plateau above r/R = 0.5. The comparison at the location x/D„ = 2.5 
indicates that the width of the predicted reverse-flow region is too thin and the predicted 
negative velocity is too small. For the turbulent shear stress shown in figure 15, the 
results of the ASE model are clearly better than those of the SKE model for all the 
locations considered. The large discrepancy seen at xj D a = 2.5 is partially due to the 
underprediction of the width of the backflow region and partially due to the experimental 
uncertainty, as evidenced by the fact that in the experimental data, the change in sign 
of the shear stress profile occurs much further away from the duct wall than the velocity 
minimum. 

4. CONCLUSIONS 

The newly developed anisotropic K-e model has been applied to calculate three conical 
diffuser flows. The comparisons with the experimental data show that the performance 
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of the new model is consistently better than that of the standard K-e model in all of the 
three flow cases. Success of the new model resides in its nonlinear stress-strain relation 
and the eddy viscosity formulation which accounts for the effect of mean shear rates. 
The new model reduces the eddy viscosity in the region of large mean strain rates where 
the standard K-e model always tends to result in overprediction. The calculations also 
indicate that the new model is as robust as the standard K-e model and requires no 
significant increase in computing time. 
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Figure 3. Turbulent kinetic energy 
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Figure 7. Displacement thickness 
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Figure 8. Momentum thickness 
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Figure 9. Shape factor 
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Figure 10. Centerline velocity 
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Figure 11. Axial mean velocity 
AKE, SKE, • Experiment 

























REPORT DOCUMENTATION PAGE 


Form Approved 
OMBNo 0704-0168 


Public reporting burden tor this collection of nformallon is estimated to average 1 hour per response, inducing the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Deectorate for Information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302. and to the Office of Management and Budget, Paperwork Reduction Project (0704-0168), Washington. DC 20503 


1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE 

October 1995 


4. TITLE AND SUBTITLE 

Calculations of Diffuser Hows With an Anisotropic K-e Model 


3. REPORT TYPE AND DATES COVERED 


Contractor Report 


S. FUNDING NUMBERS 


6. AUTHOR(S) 

J. Zhu and T.-H. Shih 


WU-505-90-5K 

NCC3-370 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESSES) 

Institute for Computational Mechanics in Propulsion 
22800 Cedar Point Road 
Cleveland, Ohio 44142 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E-9988 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESSES) 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

NASA CR-198418 

ICOMP-95-21 

CMOTT-95-4 


11. SUPPLEMENTARY NOTES 

Prepared for the 3rd International Symposium on Engineering Turbulence Modeling and Measurements, Crete, Greece, May 27-29, 
1996. J Zhu and T -H. Shih, Institute for Computational Mechanics m Propulsion and Center for Modeling of Turbulence and 
Transition, NASA Lewis Research Center ICOMP Program Director, Louis A Povmelli, orgaruzauon code 2600, (216) 433-5818 


12a. DISTRIBUTION/ AVAILABILITY STATEMENT 


12b. DISTRIBUTION CODE 


Unclassified - Unlimited 
Subject Category 34 

This publication is available from the NASA Center for Aerospace Information, (301) 621-0390 


13. ABSTRACT (Maximum 200 words) 

A newly developed anisotropic K-e model is applied to calculate three axisymmetnc diffuser flows with or without 
separation. The new model uses a quadratic stress-strain relation and satisfies the realizability conditions, i.e„ it ensures 
both the positivity of the turbulent normal stresses and the Schwarz’ inequality between any fluctuating velocities 
Calculauons arc earned out with a finite-volume method. A second-order accurate, bounded convection scheme and 
sufficiently fine gnds are used to ensure numencal credibility of the solutions. The standard K-e model is also used in 
order to highlight the performance of the new model. Companson with the experimental data shows that the anisotropic 
K-e model performs consistently better than does the standard K-e model in all of the three test cases. 


14. SUBJECT TERMS 


Turbulence model; Internal flow; Incompressible flow 


17. SECURITY CLASSIFICATION 18 SECURITY CLASSIFICATION 19 SECURITY CLASSIFICATION 
OF REPORT OF THIS PAGE OF ABSTRACT 


15 NUMBER OF PAGES 

12 


16. PRICE CODE 

A03 


20. LIMITATION OF ABSTRACT 


Unclassified 


NSN 7540-01-280-5500 


Unclassified 


Unclassified 


Standard Form 298 (Rev 2-89) 
Prescribed by ANSI Std Z39-18 
298-102 

























End of Document 



